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Abstract 

FDCSUSYDecay is a FORTRAN program package generated by FDC 
(Feynman Diagram Calculation) system fully automatically. It is ded- 
icated to calculate at tree-level all the possible 2-body decays of SUSY 
and Higgs particles in the Minimal Supersymmetric extension of the Stan- 
dard Model (MSSM). The format of its output files complies with SUSY 
Les Houches Accord and can be easily imported by other packages. 

PACS: 02.70-c; 12.60 Jv 

Program Summary 

Program Title: FDCSUSYDecay (Version 1.00) 

Journal Reference: 

Catalogue identifier: 

Keywords: SUSY decay, MSSM, FDC 

PACS: 02.70-c; 12.60 Jv 

Operating system: Linux 

Programming language: FORTRAN 77 

External libraries: CERNLIB 2003 (or up) 

Distribution format: tar gzip file 

Size of the compressed distribution file: 622,560 Bytes. 
Classification: 11.1 

Nature of physical problem: This package can calculate all the possible SUSY 
particle and Higgs 2-body decay width and branch ratio at tree-level in the 
MSSM model. 

Method of solution: By running FDC, the Feynman rules for the MSSM 
model are generated, all the decay widths are calculated analytically and corre- 
sponding FORTRAN codes are generated for this package. 



1 Generality 

With the progress in high energy physics, it is obvious that an automatic calcula- 
tion of the physical processes by the perturbative quantum field theory becomes 
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very useful and important. There are many projects for this purpose such as 
CompHEP[l], Grace[2j and so on. FDC (Feynman Diagram Calculation) [5] is 
one of such kind projects. It is mainly written in REDUCE and RLisp. As 
designed to be, it can conveniently implement physical models from the first 
principle, then construct their Lagrangians, deduce the Feynman rules and 
generate FORTRAN codes for calculations of physical quantities such as decay 
width, cross-section and matrix-element event generator. The MSSM model 
is implemented in FDC [4] and the generated results could be found on FDC 
homepag^]. To check the results such as the mixing matrices, Feynman rules, 
as well as to provide an application of FDC, we developed a program to au- 
tomatically generate the FORTRAN packages FDCSUSYDecay calculating SUSY 
particle decays and FDCEvent@LHC [5] which is a matrix-element event generator 
of the MSSM related processes at LHC. 

There are many program packages which can calculate the SUSY particle 
decays, such as multi-purpose packages PYTHIA, IS A JET and SUSYGEN or 
dedicated packages HDecay (which is special for Higgs decay in the SM and 
MSSM) [6] and SDecay [7], but the codes of these packages are manually typed. 
However, FDCSUSYDecay is a FORTRAN program package which is fully au- 
tomatically generated by FDC. In this way many uncareful errors can be sys- 
tematically debugged because these artificial errors in the model would affect 
a series of the resulting Feynman vertices and thus the corresponding decay 
equations which could be found easily. All the possible 2-body decay processes 
(in total 543 processes) for SUSY and Higgs particles in the MSSM model have 
been included. This package can calculate the decay width and branch ratio at 
Snowmass benchmark points 1-9 [8] (Tab. [J) or any user-specified parameter 
configration. The format of its output files complies with the SUSY Les Houches 
Accord (SLHA){9] and thus can be imported by other programs conveniently. 

Although presently only 2-body SUSY decays are implemented in this pack- 
age, the FORTRAN codes for any 3-body decay can be generated by FDC. Since 
the important decay channels can not be distinguished from neglectable ones 
before the MSSM parameters are chosen and to include all 3-body decays will 
make the length of the generated FORTRAN codes out of control, the 3-body 
decays are not included in the package. Alternatively, The package with needed 
3-body channels can be generated by FDC when the MSSM parameters fixed. 

The package can be used in two modes. The first mode is called the high- 
energy scale mode. In this mode, users need to input the parameters in high- 
energy scale scenarios (mSUGRA, GMSB and AMSB) or choose from nine Snow- 
mass benchmark points (which have been fed into this package hard-wiredly). 
Then the program obtains all the parameters at SUSY-breaking scale through 
interface with IS A JET [10] and then calculates everything such as masses, mix- 
ing matrices, coupling constants and decay width. The second mode is called 
low-energy scale mode, in which users need to input manually the MSSM pa- 
rameters: 

• Mi,M 2 ,M 3 : U(l), U(2), U(3) gaugino mass parameters 
^httpTTTww . ihep . ac . cn/limuen/uj x/public_html 
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Table 1: High-energy scale scenario points corresponding to Snowmass points 

i-9[n] 



• : Higgsino mass parameter 

• tan /3 : the ratio of the vev's of the two-Higgs doublet fields 

• A T ,Ai,,At : the 3rd generation trilinear couplings 

• m A o : the neutral pseudoscalar Higgs (^4°) mass 

• m? (= m 2 ),m? , m 2 (= m? ),m? , 

m l (= m 2 - ),m~ , m 2 ~ , m 2 ~ (= m? ),m~ , m? : the SUSY-breaking mass param- 
eters for the 1st /2nd generation sfermions 



'xi ,m? , m 2 . f— ™ 2 
t 2 ' b 2 ' ii v 

generation sfermions 



• mi ,mi ,m- ,mi (= mi ),mi : the SUSY-breaking mass parameters for the 3rd 

n t 2 bo t 1 v bi t 2 ° 



The Yukawa couplings for the 1st and 2nd generation fermions are fixed to zero 
as an approximation and the corresponding soft SUSY breaking terms are also 
fixed to zero in the MSSM model used here. If needed, they can be easily 
included in our package by FDC. 



2 Installation and Usage 
2.1 Installation 

The package is distributed in a gzip compressed file named fsdecay.tar.gz and 
can be downloaded from FDC homepage: 

http : //www. ihep . ac . cn/lunwen/wjx/public_html 
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The CERN program library (CERNLIB) 2003 Ql] (or up) and ISAJET are 
required only for high-energy scale mode. The installation steps are given below 
and for convenience it is assumed 

1. CERNLIB is located at directory /cern/2003, 

2. the current user account is linuxuser; 

3. the downloaded file fsdecay.tar.gz are stored at directory ~linuxuser; 

4. the package will be installed to ~linuxuser/FDCDecay. 
The installation steps are 

1. Create the target directory and go into it by 

mkdir "linuxuser/FDCDecay 
cd "linuxuser/FDCDecay 

2. Decompress the downloaded file by 

tar -xzvf ~linuxuser/f sdecay .tar .gz 

3. Specify the CERNLIB path by adding following line to the very begin- 
ning of the file makefil^ 

cern = /cern/2003 

4. Compile the low-energy scale mode by 

make lowscale 

5. Compile the high-energy scale mode by 

make highscale 

6. Delete the temporary files by 

make clean 

After steps above without encountering any errors, one will obtain two exe- 
cutable files named fsdecayhigh and fsdecaylow in directory ~linuxuser/FDCDecay 
which are stand-alone programs used to calculate decay width and branch ratio 
in high-energy scale and low-energy scale modes respectively. 

2 Users can substitute the user account and paths for their owns, when installing on their 
own computers. 

3 This step is necessary only for high-energy scale mode. 
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2.2 Usage 



As mentioned above, with this package, users can calculate decays in high-energy 
scale mode or low-energy scale mode. While in high-energy scale mode, users 
can either specify the scenarios (including mSUGRA, GMSB, AMSB) and their 
parameters manually or just use the Snowmass benchmark points. While in 
low-energy scale mode, users must input the MSSM (and the Standard Model) 
parameters manually. 

For this purpose, we prepared two executable programs fsdecayhigh and 
fsdecaylow corresponding to these two modes respectively. 

Next we will describe the usage of these two programs. Test-run examples 
are given in Sec. [6l 

2.2.1 Using high-energy scale mode 

It is very simple to calculate decay width in high-energy scale mode. What 
users need to do is just run program fsdecayhigh and then the program will give 
prompts on screen and require corresponding inputs. 

2.2.2 Using low-energy scale mode 

To calculate decay width in this mode is somewhat more complicated. Users 
should follow steps below: 

1. Modify the FORTRAN file usrpara.f in the installation directory to assign 
value to parameters. There is a subroutine FUPAR in which users can 
specify the MSSM and SM parameters respectively by assigning values to 
corresponding variables. All these variables are listed in TabfSJ 

2. Re-build the executable program fsdecaylow by command: 
make lowscale clean 

3. Run program fsdecaylow and it will give prompt for users to specify the 
output file name 



3 Output files 

The format of the generated output file for decays complies with the SUSY Les 
Houches Accord|9j. It is composed of several blocks which are separated by 
empty lines. Eacn of these blocks describes a SUSY particle decay. Every block 
begins with a DECAY statement specifying the decaying mother particle and its 
total width, in the format as follows 

# PDG WIDTH(GeV) 

DECAY 2000006 9 . 00247665E+00 # ~t_2 decays 
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Table 2: Variables corresponding to the MSSM and SM parameters in subrou- 
tine FUPAR 



The first integer is the PDG particle code[13] of the decaying particle. The 
subsequent real number is its total width. The end comment gives a human 
readable translation of the PDG code. All the decay modes for this mother 
particle are listed in the subsequent lines in the format 
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In each line above, the first real number is the branch ratio of the decay 
mode, the first integer is the number of daughters (which is always 2 in this 
package) and the following integers are the PDG codes of these daughters. The 
specific FORTRAN format for the DECAY statement and the entries in a block 
are respectively 

FORMAT ( ' # ' , 7X , ' PDG ' , 9X , 'WIDTH(GeV) ' ) 
FORMAT ( 'DECAY' , IX , 19 ,3X , IP ,E16 . 8 , OP , 3X , '# ' , IX, A25) 
FORMAT ( '# ' ,4X , 'BR' , 16X, ' NDA ' , 3X , ' ID1 ' , 7X , ' ID2 ' ) 
F0RMAT(3X,1P,E16.8,0P,3X,I2,3X,2(I9,1X) ,2X, '#' ,1X,A25) 



4 Cross-Check with ISAJET 

4.1 Running mass parameters and higher-order mass cor- 
rections 

In ISAJET (Version 7.64), running mass for t, b and t are used in the calculations 
of stau(fi), sbottom (bi) and stop (ti) mass eigenvalues, mixing angles and in 
the vertices of j ermionir jb j t) ~ sfermion(fi/bi/ti) ~ gaugino(xf / xl) ■ 
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Table 3: Higgs and gluino mass in ISAJET (I) and FDCSUSYDecay (F) 



Higher order corrections are added in ISAJET to the mass of Higgs H°, h° , 
H + and gluino g. The corrections at SPS 1-9 bring less than 1% modification 
to the H and H + mass, 27% ~ 43% to h° mass and 5% ~ 15% for g mass. 
In Tabf3]we listed the Higgs and gluino mass at SPS 1-9 from ISAJET and our 
decay package. However these corrections have not been implemented in FDC 
yet (thus in our package). 

In order to compare the decay widths, we take, temporarily, the running 
mass and the corrected Higgs and gluino mass as inputs from ISAJET. 

4.2 Cross-check and results 

Taking higher order corrections to Higgs and gluino mass as inputs, we cal- 
culated the decay widths at SPS 1-9 and compared them with the data from 
ISAJET. 

For the decay modes with 10% discrepancies or larger between ISAJET and 
our package, we analyzed the decay equations in the source codes of ISAJET 
(Version 7.64). Two reasons for the discrepancies were found: (1) correction to 
triple Higgs vertices and QCD radiative correction to Higgs — > q(t/b) + q(t/b) 
included in ISAJET; (2) different decay equations used by ISAJET and our 
package which imply something wrong in the codes of ISAJET0. 

For the first kind of discrepancy, we list the widths of Higgs decays which 
are influenced by the higher order corrections in Tab[4l From Tab|4]we can see 
that the corrections bring large modifications sometimes. 

In the following we present a detailed disscussion for the second kind of 
discrepancy. There are 3 kinds of decay modes with different decay equations 
between ISAJET and FDCSUSYDecay. They are: 

1. A chargino decays into a left-sstrange and a c-quark (xf — > sl + c); 

2. The neutral pseudoscalar Higgs decays into 2 different staus (A — > f-f + 

4 We carefully checked the MSSM model used in FDC with that in Ref. [U] 
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Table 4: Higgs decay widths in ISAJET (I) and FDCSUSYDecay (F) 
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Table 5: Decay width from ISAJET (I) and FDCSUSYDecay (F) 



3. The neutral pseudoscalar Higgs decays into a chargino pair (A — > \i + 
X-.i = h'2). 

The numerical results of decay widths at SPS 1 to 9 are shown in Tab[5]and it 
is very clear that there are large discrepancies in these decay channels. 

In the following, we briefly present the definition of the particles, the result- 
ing mixing matrices and the related Feynman rules. The mass matrix for stau 
in the tl — tr basis is 



M f = 
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m T (A T - //tan/3) M$ 



2 



M? L = m\ + m T1 + m 2 cos2/3(-i + sin 2 6 W ), (1) 

M~ R = m 2 + fh\ - m 2 z cos 2/3 sin 2 8 W , (2) 

where m~ , m? 2 and A T are parameters in soft SUSY breaking term. The stau 
mass matrix is diagonalized by a unitary matrix U T 
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We adopt the positively charged charginos as particles and define the physical 



states with mass m^ as 



*(\/ )= ( ^ ), *(x7) = *(xtf= ( A | j {->) 
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where 



The A ± denote winos and Hy and denote charged higgsinos. The unitary 
matrices U and V are given as 



sin cos9p J \ r]2 J ' V — sin6>£ cos 0£ 

The chargino mass matrix is 

M _ = f _ M 2 V2m z cos9 w smf3 \ 

x \ y/2m z cos #w cos /3 /i / ' 

It can be diagonalized by unitary matrices in ^ 

U-'M.V = ( ^ ^ ) , (|m £2 | > |m gl |), (8) 

where rji =sign(m£ 1 ) and r]2 =sign(mg 2 ). The chargino mass then will be 
m xl =Vi m di, m x2 = mm d2 . (9) 

The related Feynman rules used in our package are given in FigQ] and El 
which have been carefully checked with those in Ref. [14]. 
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Figure 1: Vertices of X^~Sz,c in FDC 
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and decay width equations are listed as: 
mt^sz+c) = ^^(m 2 -<+ m|J cos 2 fl H (10) 
r(x^^ + c) = ^^(m 2 -m 2 [+ m 2 2 )sin 2 ^ (11) 
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r(^ff+f+) = ° m i Pa i (H + A T tan/3) 2 
327rm^oTO^ 

2|p I 

r(A°^x^+Xi") = ' cm| (cos/3cos6> fl sin6> L +sin/3sin6' fl cos6'z,) 2 13) 

2| p I 

r(^ ^X2"+X^) = ' (cos /3 sin 6 R cos 6> L + sin cos 6> fl sin 0^)^14) 

where |P cm | is the momentum of one of the final particle in the rest frame of 
the initial particle 

_ y/[Af 2 - (mi + m 2 ) 2 ] [Af 2 - (mi - m 2 ) 2 ]' 
1 cm| ~ 2M 

where M is the mass of the initial particle and m\,m% is the masses of the final 
particles. 
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5 Summary 



We have presented the FORTRAN program package FDCSUSYDecay, which 
is used to calculate all the possible SUSY particle and Higgs 2-body decays 
in the MSSM model at tree-level. The package with needed 3-body channels 
can be generated by FDC when the MSSM parameters are fixed. One of the 
attractive features of this package is that all the FORTRAN codes are generated 
by FDC with full automation and it can avoid uncareful errors. Although this 
package can now only work at tree-level, the ability of one-loop calculation in 
FDC are under developing. A new version of FDCSUSYDecay with full one-loop 
correction will be generated in near future. 



6 Test Run and Output 

6.1 Test-run in high-energy scale mode with parameters 
of Snowmass benchmark point 1 

The screen output when running fsdecayhigh is as following: 

* FDCSUSYDecay 1.00 * 
Enter Snoumass point (1-9): 

OR Enter to specify scenario and parameters manually: 
1 

Enter output file name (in single quotes) : 
'testrunl . dat ' 



The decay table is stored in file testrunl.dat and reads like (partly listed): 

# ===================================== 

# DECAY TABLE CREATED BY FDCSUSYDECAY 

# DATE: 11/ 4/2006 16:43:46 

# ===================================== 
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~t. 


.2 


-> 


t 


"chiO. 


.1 



# PDG WIDTH(GeV) 

DECAY 1000021 6 . 33793298E+00 # ~g decays 

# BR NDA ID1 ID2 
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1 , 


, 1922412ohi-01 


2 


5 


-1 A AA AA C 

-100000b 


# 


~g 


- > 


b "*b_ 


.lbar 


1 . 


, 1922412oh-01 


2 


-5 


-1 A AA AA C 

1000005 


# 


~g 


-> 


bar 


b_l 


5 . 


, 10bbllooh-02 


2 


2 


AA AA AA A 

-2000002 


# 


~g 


- > 


u ~u_ 


_Rbar 


5 . 


, 10bbllooE-02 


2 


-2 


AA AA AA A 

2000002 


# 


"g 


-> 


ubar 


~u_R 


5 . 


H AC /I A AC "VC A A 

, 10b429b fh-02 


2 


4 


AA AA AA / 

-2000004 


# 


~g 


- > 


c ~c_ 


.Rbar 


5 , 


H AC /I A AC "VC A A 

, 10b429b fh-02 


2 


-4 


AA AA AA / 

2000004 


# 


"g 


- > 


cbar 


"c_R 


5 . 


AOAjCATCAT -1 aa 

. 08960759E-02 


2 


1 


AA AA AA -1 

-2000001 


# 


"g 


-> 


d ~d_ 


_Rbar 


5 , 


A O A£? ATC AC A A 

, 0o9b0 ( b9E-02 


2 


- 1 


AA AA AA -1 

2000001 


# 


~g 


- > 


dbar 


"d_R 


g 


OfiQRfil Q7F 09 


2 


3 






g 


_ > 




Rbar 


5, 


08958197E-02 


2 


-3 


2000003 


# 


~g 


-> 


sbar 


~s_R 


3. 


90571495E-02 


2 


6 


-1000006 


# 


"g 


-> 


t "t. 


.lbar 


3. 


.90571495E-02 


2 


-6 


1000006 


# 


~g 


-> 


tbar 


~t_l 


3, 


89890226E-02 


2 


5 


-2000005 


# 


~g 


-> 


b ~b. 


.2bar 


3 


.89890226E-02 


2 


-5 


2000005 


# 


"g 


-> 


bar 


"b_2 


2. 


81103874E-02 


2 


2 


-1000002 


# 


~g 


-> 


u ~u. 


.Lbar 


2. 


.81103874E-02 


2 


-2 


1000002 


# 


"g 


-> 


ubar 


~u_L 


2. 


■80978313E-02 


2 


4 


-1000004 


# 


"g 


-> 


c "c. 


.Lbar 


2. 


80978313E-02 


2 


-4 


1000004 


# 


"g 


-> 


cbar 


~c_L 


2, 


, 13047293E-02 


2 


1 


-1000001 


# 


~g 


-> 


d ~d. 


.Lbar 


2. 


. 13047293E-02 


2 


-1 


1000001 


# 


"g 


-> 


dbar 


~d_L 


2. 


. 13044513E-02 


2 


3 


-1000003 


# 


"g 


-> 


s "s. 


.Lbar 


2, 


, 13044513E-02 


2 


-3 


1000003 


# 


~g 


-> 


sbar 


~s_L 



# PDG WIDTH(GeV) 

DECAY 1000003 5 . 00062352E+00 # ~s_L decays 



BR 


NDA 


ID1 


ID2 












6.10919249E-01 


2 


4 


-1000024 


# 


~s_L 


-> c 


~chi-_ 


.1 


3.10762902E-01 


2 


3 


1000023 


# 


~s_L 


-> s 


"chiO. 


.2 


3.89466555E-02 


2 


4 


-1000037 


# 


~s_L 


-> c 


~chi-. 


.2 


2.41598639E-02 


2 


3 


1000022 


# 


~s_L 


-> s 


"chiO. 


.1 


1.36997584E-02 


2 


3 


1000035 


# 


~s_L 


-> s 


"chiO. 


.4 


1.51157142E-03 


2 


3 


1000025 


# 


~s_L 


-> s 


"chiO. 


.3 



6.2 Test-run in high-energy scale mode with manually in- 
puted parameters 

The screen output when running fsdecayhigh is as following: 

* FDCSUSYDecay 1.00 * 
Enter Snoumass point (1-9): 

OR Enter to specify scenario and parameters manually: 


1. mSUGRA 

2. minimal Gauge-Mediated SUSY Breaking (mGMSB) 

3. Anomaly-Mediated SUSY Breaking (AMSB) 
Please choose scenario (1-3) : 

1 

Enter M_0, M_(l/2), A_0, tan(BETA) , sign(MU) , M_t(top mass) 
100.0 260.0 -100.0 12.0 1.0 175.0 
Enter output file name (in single quotes): 
'testrun2 . dat ' 
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The decay table is stored in file testrun2.dat and reads like (partly listed): 



# DECAY TABLE CREATED BY FDCSUSYDECAY 

# DATE: 11/ 4/2006 16:50:38 

# ===================================== 

# 

# PDG WIDTH(GeV) 





1000021 


6 


.87489451E+00 




ft ~g decays 












BR 




NDA 


ID1 




ID2 












1, 


22817522E- 


■01 


2 




5 


-1000005 


# 


O 


-> 


b ~b. 


.lbar 


1. 


22817522E- 


■01 


2 




-5 


1000005 


# 


"g 


-> 


bar 


'b_l 


5, 


■00789610E- 


■02 


2 




1 


-2000001 


# 


~g 


-> 


d ~d. 


.Rbar 


5. 


.00789610E- 


-02 


2 




-1 


2000001 


# 


"g 


-> 


dbar 


~d_R 


5, 


,00787337E- 


■02 


2 




3 


-2000003 


# 


~g 


-> 


s ~s_ 


.Rbar 


5. 


■00787337E- 


-02 


2 




-3 


2000003 


# 


"g 


-> 


sbar 


~s_R 


4, 


98299423E- 


■02 


2 




2 


-2000002 


# 


~g 


-> 


u ~u_ 


.Rbar 


4. 


98299423E- 


■02 


2 




-2 


2000002 


# 


~g 


-> 


ubar 


~u_R 


4. 


.98194437E- 


-02 


2 




4 


-2000004 


ft 


"g 


-> 


c "c. 


.Rbar 


4. 


.98194437E- 


■02 


2 




-4 


2000004 


# 


"g 


-> 


cbar 


~c_R 


4. 


34327716E- 


-02 


2 




6 


-1000006 


# 


"g 


-> 


t ~t_ 


.lbar 


4, 


34327716E- 


■02 


2 




-6 


1000006 


# 


~g 


-> 


tbar 


~t_l 


3. 


■66751775E- 


-02 


2 




5 


-2000005 


# 


"g 


-> 


b "b. 


.2bar 


3, 


.66751775E- 


■02 


2 




-5 


2000005 


# 


"g 


-> 


bar 


*b_2 


2. 


74132411E- 


■02 


2 




2 


-1000002 


# 


~g 


-> 


u ~u. 


.Lbar 


2, 


74132411E- 


-02 


2 




-2 


1000002 


ft 


"g 


-> 


ubar 


~u_L 


2. 


,74020784E- 


■02 


2 




4 


-1000004 


# 


"g 


-> 


c "c. 


.Lbar 


2. 


,74020784E- 


-02 


2 




-4 


1000004 


# 


"g 


-> 


cbar 


~c_L 


2, 


, 12261878E- 


■02 


2 




1 


-1000001 


# 


"g 


-> 


d ~d. 


.Lbar 


2. 


, 12261878E- 


-02 


2 




-1 


1000001 


# 


"g 


-> 


dbar 


~d_L 


2. 


. 12259409E- 


■02 


2 




3 


-1000003 


# 


"g 


-> 


s ~s_ 


.Lbar 


2. 


. 12259409E- 


-02 


2 




-3 


1000003 


# 


"g 


-> 


sbar 


~s_L 



# PDG WIDTH(GeV) 



DECAY 2000006 


9.48857394E+00 # 


~t_2 decays 














# BR 


NDA 


ID1 


ID2 
















2.19460985E-01 


2 


5 


1000037 


# 


"t. 


.2 


-> 


b 


~chi+. 


.2 


2.12621918E-01 


2 


23 


1000006 


# 


~t. 


.2 


-> 


ZO 


"t_l 




1.97473598E-01 


2 


6 


1000035 


# 


"t. 


.2 


-> 


t 


"chiO. 


.4 


1.65481446E-01 


2 


5 


1000024 


# 


"t. 


.2 


-> 


b 


~chi+. 


.1 


6.29074190E-02 


2 


6 


1000023 


# 


"t. 


.2 


-> 


t 


"chiO. 


.2 


4.22399775E-02 


2 


6 


1000025 


# 


"t. 


.2 


-> 


t 


"chiO. 


.3 


4.13479750E-02 


2 


24 


1000005 


# 


~t. 


.2 


-> 


W+ 


"b_l 




3.30474133E-02 


2 


25 


1000006 


# 


"t. 


.2 


-> 


hO 


~t_l 




2.54192679E-02 


2 


6 


1000022 


# 


~t. 


.2 


-> 


t 


"chiO. 


.1 



# 

DECAY 



PDG 

1000003 
BR 

6.11686309E-01 
3.11939618E-01 
3.86054484E-02 
2.27239525E-02 
1.35532277E-02 
1.49144465E-03 



WIDTH(GeV) 
5.17658881E+00 

NDA ID1 

2 

2 

2 

2 

2 

2 



~s_L decays 

ID2 
-1000024 

1000023 
-1000037 

1000022 

1000035 

1000025 



# - 


s_L 


-> c 


~chi-_ 


.1 


# - 


s_L 


-> s 


"chiO. 


.2 


# - 


s_L 


-> c 


"chi-. 


.2 


# - 


s_L 


-> s 


"chiO. 


.1 


# - 


s_L 


-> s 


"chiO. 


.4 


# - 


s_L 


-> s 


"chiO. 


.3 
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6.3 Test-run in low-energy scale mode 

The file usrpara.f is prepared as: 

SUBROUTINE FUPARO 

c 

c By this subroutine, users can specify the MSSM and SM 

c parameters at low-energy scale. 

c 

c ** IMPORTANT ** 

c PROGRAM DO NOT CHECK THE VALIDITY OF PARAMETERS VALUE 

c YOU INPUTED. 

c 

c **REMEMBER TO RE-BUILD THE PROGRAM AFTER EVERY MODIFICATION 

c BY COMMAND: 

c make louscale 
c 

c 

c The MSSM parameter variables are listed below: 

c The unit is GeV 

c 

c FSSGM1 = U(l) gaugino mass parameter 

c FSSGM2 = U(2) gaugino mass parameter 

c FSSGM3 = U(3) gaugino mass parameter 

c FSSHMU = Higgsino mass parameter 

c FSSTANB = ratio of Higgs vacuum expectation 

c FSSMHA = neutral pseudoscalar Higgs (AO) mass 

c FSSAAL = trilinear coupling for tau lepton 

c FSSAAT = trilinear coupling for t-quark 

c FSSAAB = trilinear coupling for b-quark 

c FSSSBEL = SUSY-breaking mass parameter for left s-electron 

c FSSSBER = SUSY-breaking mass parameter for right s-electron 

c FSSSBMUL = SUSY-breaking mass parameter for left s-muon 

c FSSSBMUR = SUSY-breaking mass parameter for right s-muon 

c FSSSBTAU1 = SUSY-breaking mass parameter for left s-tau 

c FSSSBTAU2 = SUSY-breaking mass parameter for right s-tau 

c FSSSBUPL = SUSY-breaking mass parameter for left s-up 

c FSSSBUPR = SUSY-breaking mass parameter for right s-up 

c FSSSBCHL = SUSY-breaking mass parameter for left s-charm 

c FSSSBCHR = SUSY-breaking mass parameter for right s-charm 

c FSSSBDNR = SUSY-breaking mass parameter for right s-down 

c FSSSBSTR = SUSY-breaking mass parameter for right s-strange 

c FSSSBBT2 = SUSY-breaking mass parameter for heavy s-bottom 

c FSSSBTP1 = SUSY-breaking mass parameter for light s-top 

c FSSSBTP2 = SUSY-breaking mass parameter for heavy s-top 

c 

c The SM parameter variables are listed below: 

c The unit is GeV 

c 

c FSMME = electron mass 

c FSMMMU = muon mass 

c FSMMTAU = tau mass 

c FSMMUP = u-quark mass 

c FSMMDN = d-quark mass 

c FSMMCH = c-quark mass 
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FSMMST 
FSMMTP 
FSMMBT 
FSMMZ 



s-quark mass 
t-quark mass 
b-quark mass 
Z boson mass 



FSMSN2THW = ueak-mixing angle SIN~2(THETA_W) 
FSMFNST = fine-structure constant ALPHA_EM 
FSMALPS = strong-coupling constant ALPHA_S 



IMPLICIT NONE 

c BLOCK FUSSPAR START 

DOUBLE PRECISION FSSGM1 , FSSGM2 , FSSGM3 , FSSHMU , FSSTANB , FSSMHA , 
$ FSSAAL , FSSAAT , FSSAAB , FSSSBEL , FSSSBER , FSSSBMUL , FSSSBMUR , 

$ FSSSBTAU1,FSSSBTAU2,FSSSBUPL,FSSSBUPR,FSSSBCHL,FSSSBCHR, 
$ FSSSBDNR,FSSSBSTR,FSSSBBT2,FSSSBTP1,FSSSBTP2 

COMMON /FUSSPAR/ FSSGM1 , FSSGM2 ,FSSGM3 , FSSHMU, FSSTANB, FSSMHA, 
$ FSSAAL , FSSAAT , FSSAAB , FSSSBEL , FSSSBER , FSSSBMUL , FSSSBMUR , 

$ FSSSBTAU1,FSSSBTAU2,FSSSBUPL,FSSSBUPR,FSSSBCHL,FSSSBCHR, 
$ FSSSBDNR,FSSSBSTR,FSSSBBT2,FSSSBTP1,FSSSBTP2 
c BLOCK END 

c BLOCK FUSMPAR START 

DOUBLE PRECISION FSMME , FSMMMU , FSMMTAU , FSMMUP , FSMMDN , FSMMCH , 
$ FSMMST , FSMMBT , FSMMTP , FSMMZ , FSMSN2THW , FSMFNST , FSMALPS 

COMMON /FUSMPAR/ FSMME, FSMMMU, FSMMTAU, FSMMUP, FSMMDN, FSMMCH, 
$ FSMMST , FSMMBT , FSMMTP , FSMMZ , FSMSN2THW , FSMFNST , FSMALPS 
c BLOCK END 

c SET YOUR MSSM PARAMETERS HERE 



FSSGM1 




99.6476898d0 


FSSGM2 




192.514954d0 


FSSGM3 




579.09314d0 


FSSHMU 




354.201385d0 


FSSTANB 




10.0DO 


FSSMHA 




394.874817D0 


FSSAAL 




-253.15593d0 


FSSAAB 




-766.899109d0 


FSSAAT 




-495.786896d0 


FSSSBEL 




38581. 48830D0 


FSSSBER 




18547. 65040D0 


FSSSBMUL 




38581. 48830D0 


FSSSBMUR 




18547. 65040D0 


FSSSBTAU1 




38243. 67190D0 


FSSSBTAU2 




17847. 84380D0 


FSSSBUPL 




289769. 750D0 


FSSSBUPR 




271373. 5620D0 


FSSSBCHL 




289769. 750D0 


FSSSBCHR 




271373. 5620D0 


FSSSBDNR 




269598. 8440D0 


FSSSBSTR 




269598. 8440D0 


FSSSBBT2 




266551. 1560D0 


FSSSBTP1 




243069. 0310D0 


FSSSBTP2 




178026. 6410D0 
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cc SET YOUR SM PARAMETERS HERE 



FSMME 




5.10999991D-4 


FSMMMU 




1.04999997D-1 


FSMMTAU 




1.77699995 


FSMMUP 




5.5999998D-3 


FSMMDN 




9.89999995D-3 


FSMMCH 




1.35000002D0 


FSMMST 




1.99000001D-1 


FSMMTP 




175D0 


FSMMBT 




5D0 


FSMMZ 




91.1699982D0 


FSMSN2THW 




0.2311678d0 


FSMFNST 




1D0/127. 970535 


FSMALPS 




0.118170463d0 



END 

The screen output when running fsdecaylow is as following: 

# FDCSUSYDecay 1.00 * 

ENTER output file name (in single quotes): 
'testrun3 . dat ' 

The decay table is stored in file testrun3.dat and reads like (partly listed): 

# ====================================== 

# DECAY TABLE CREATED BY FDCSUSYDECAY 

# DATE: 11/ 4/2006 16:58:46 

# ====================================== 

### 

# PDG WIDTH(GeV) 



DECAY 2000006 < 


K00247595E+00 # 


"t_2 dec 


ays 














# BR 


NDA 


ID1 


ID2 
















2.25878090E-01 


2 


23 


1000006 


# 


~t. 


.2 


-> 


ZO 


"t_l 




2.21288384E-01 


2 


5 


1000037 


# 


~t. 


.2 


-> 


b 


~chi+. 


.2 


1.90569924E-01 


2 


6 


1000035 


# 


"t. 


.2 


-> 


t 


"chiO. 


.4 


1.61424752E-01 


2 


5 


1000024 


# 


"t. 


.2 


-> 


b 


~chi+. 


.1 


6.33183886E-02 


2 


6 


1000023 


# 


"t. 


.2 


-> 


t 


"chiO. 


.2 


3.98983144E-02 


2 


6 


1000025 


# 


"t. 


.2 


-> 


t 


"chiO. 


.3 


3.91332354E-02 


2 


24 


1000005 


# 


~t. 


.2 


-> 


W+ 


"b_l 




3.31751871E-02 


2 


25 


1000006 


# 


~t. 


.2 


-> 


hO 


"t_l 




2.53137244E-02 


2 


6 


1000022 


# 


"t. 


.2 


-> 


t 


"chiO. 


.1 



# PDG WIDTH(GeV) 





1000021 


6 


.33793301E+00 


# ~g decays 












BR 




NDA 


ID1 


ID2 












1. 


, 19224123E- 


-01 


2 


5 


-1000005 


# 


~g 


-> 


b "b. 


.lbar 


1, 


, 19224123E- 


■01 


2 


-5 


1000005 


# 


~g 


-> 


bar 


"b_l 


5. 


. 10661146E- 


-02 


2 


2 


-2000002 


# 


"g 


-> 


u "u. 


.Rbar 


5, 


, 10661146E- 


-02 


2 


-2 


2000002 


# 


~g 


-> 


ubar 


~u_R 


5. 


. 10542966E- 


-02 


2 


4 


-2000004 


# 


"g 


-> 


c "c. 


.Rbar 


5. 


, 10542966E- 


-02 


2 


-4 


2000004 


# 


~g 


-> 


cbar 


~c_R 


5. 


■08960753E- 


-02 


2 


1 


-2000001 


# 


"g 


-> 


d ~d. 


.Rbar 


5. 


08960753E- 


-02 


2 


-1 


2000001 


# 


"g 


-> 


dbar 


~d_R 


5. 


.08958191E- 


-02 


2 


3 


-2000003 


# 


~g 


-> 


s ~s_ 


.Rbar 


5. 


■08958191E- 


-02 


2 


-3 


2000003 


# 


"g 


-> 


sbar 


"s_R 
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3. 


.90571491E-02 


2 


6 


-1000006 


# 


~g 


-> 


t "t. 


.lbar 


3. 


.90571491E-02 


2 


-6 


1000006 


# 


~g 


-> 


tbar 


"t_l 


3. 


89890229E-02 


2 


5 


-2000005 


# 


~g 


-> 


b ~b. 


_2bar 


3. 


89890229E-02 


2 


-5 


2000005 


# 


"g 


-> 


bar 


"b_2 


2, 


.81103873E-02 


2 


2 


-1000002 


# 


~g 


-> 


u ~u. 


.Lbar 


2. 


81103873E-02 


2 


-2 


1000002 


# 


~g 


-> 


ubar 


~u_L 


2. 


.80978312E-02 


2 


4 


-1000004 


# 


"g 


-> 


c "c. 


_Lbar 


2, 


80978312E-02 


2 


-4 


1000004 


# 


~g 


-> 


cbar 


~c_L 


2. 


13047294E-02 


2 


1 


-1000001 


# 


"g 


-> 


d ~d. 


_Lbar 


2. 


, 13047294E-02 


2 


-1 


1000001 


# 


~g 


-> 


dbar 


"d_L 


2. 


, 13044514E-02 


2 


3 


-1000003 


# 


"g 


-> 


s ~s. 


_Lbar 


2. 


, 13044514E-02 


2 


-3 


1000003 


# 


~g 


-> 


sbar 


~s_L 



# PDG WIDTH(GeV) 

DECAY 1000003 5 . 00062300E+00 # ~s_L decays 



BR 


NDA 


ID1 


ID2 














6.10919249E-01 


2 


4 


-1000024 


# 


~s_L 


-> 


c 


~chi-_ 


.1 


3.10762900E-01 


2 


3 


1000023 


# 


~s_L 


-> 


s 


"chiO. 


.2 


3.89466552E-02 


2 


4 


-1000037 


# 


~s_L 


-> 


c 


~chi-_ 


.2 


2.41598663E-02 


2 


3 


1000022 


# 


~s_L 


-> 


s 


"chiO. 


.1 


1.36997584E-02 


2 


3 


1000035 


# 


~s_L 


-> 


s 


"chiO. 


.4 


1.51157143E-03 


2 


3 


1000025 


# 


~s_L 


-> 


s 


"chiO. 


.3 
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